Gas-Mediated Impact Dynamics in Fine-Grained Granular Materials 
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Non-cohesive granular media exhibit complex responses to sudden impact that often differ from 
those of ordinary solids and liquids. We investigate how this response is mediated by the presence 
of interstitial gas between the grains. Using high-speed x-ray radiography we track the motion of 
a steel sphere through the interior of a bed of fine, loose granular material. We find a crossover 
from nearly incompressible, fluid-like behavior at atmospheric pressure to a highly compressible, 
dissipative response once most of the gas is evacuated. We discuss these results in light of recent 
proposals for the drag force in granular media. 

PACS numbers: 45.70.Cc, 47.56.+r, 83.80.Fg, 83.10.Tv 



Studies of the impact of solid objects into granular beds 
have a long history, with first systematic work dating 
back to the 18th century [1, 2]. The topic reemerged in 
the 1960's with attempts to characterize the strength of 
the lunar surface |3|, [Jl and, recently, with investigations 
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161 , im . One unresolved issue remains 
the form for the drag force experienced by the impacting 
solid as it moves through the granular medium. Vari- 
ous competing force laws have been proposed, based on 
scaling laws relating the penetration depth to the impact 
energy or momentum [s', o], or on direct measurement of 
the trajectory of the impacting object 0, H, 0, El- 
However, none of these take into account the interaction 
between the solid grains inside a granular bed and the 
surrounding gas. 

This interaction is known to play a significant role in 
many situations where the bed is externally forced, espe- 
cially for small grain sizes, when the bed's gas permeabil- 
ity becomes sufficiently low to sustain a pressure gradient 
that can compete with the weight of the material. The 
resulting feedback between grain motion and ambient gas 
flow gives rise to complex dynamics not only when a bed 
is fluidized by direct gas injection Il8^, but also in many 
vibrated granular systems [1^, HO, Ull, |22[. There have 
been indications that the morphology of craters differs 
for small grain sizes due to interstitial gas flows |13( . 
Also, for fine powders, the size and shape of granular 
jets ejected upward after impact was found to depend on 
the presence of interstitial gas flG*]. However, with one 
exception [17|] , the role of gas-grain interactions in deter- 
mining the trajectory of impacting solids has not been 
investigated in detail. Here we demonstrate that ambi- 
ent gas inside a granular bed strongly affects the impact 
dynamics and show how this alters the drag force in a 
fine-grained granular medium. 

Previous experiments used either two-dimensional 
(2D) set-ups [Uj or indirect methods to gauge the motion 
of a projectile inside a 3D bed [I, @, 0, [1, @, 11Q|. In con- 
trast, our approach is based on high-speed x-ray imaging 
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FIG. 1: (Color Online). Pressure dependence of impact dy- 
namics. Composite x-ray images at 12 kPa (a) 5 ms, (b) 40 
ms and (c) 57 ms after impact. Images at 0.7 kPa (d) 5 ms, 
(e) 12 ms and (f) 28 ms after impact. 



which gives direct, time-resolved access to the dynam- 
ics in the bed interior and also allows us to extract local 
changes in the bed packing density. 

For the experiments reported here, a steel sphere (di- 
ameter Dg = 12mm) was dropped from a height of 0.34 
m into an 85 mm deep bed of boron carbide (B4C) parti- 
cles (50 /im avg. diameter). B4C, which is non-spherical, 
was chosen to optimize the x-ray transmission; separate 
experiments, studying jet formation in a variety of dif- 
ferent media, indicate that grain shape is not a critical 
parameter 1^, 2^- The bed was contained in a cylindri- 
cal tube with 35 mm inner diameter. Before each drop 
the bed was aerated by dry nitrogen entering through 
a diffuser built into the bottom of the container. By 
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slowly turning off the nitrogen flow, the packing fraction 
(j) = Vg/Vtot^ where Vg is volume occupied by grains and 
Vtot is the total volume of the bed, was adjusted before 
each drop to a value around 0.5. The system could be 
sealed and evacuated down to pressures as low as 0.7 
kPa. The pump speed was limited to prevent air flow 
from disturbing the loose packing. We checked for elec- 
trostatic charging by performing experiments in air at 
a high humidity (~50%) where electrostatic effects typ- 
ically vanish [24] and observed no qualitative change in 
the impact dynamics. 

X-ray imaging was performed at the GSECARS beam- 
line at the Advanced Photon Source using a high inten- 
sity beam with energy width 5 keV centered at 22 keV. 
X-ray transmission through the bed was imaged off a 
phosphor screen at 6000 frames per second using a Phan- 
tom v7 video camera. The beam size restricted the fleld 
of view to 22 mm x 8.7 mm sections of the bed. To 
capture the dynamics across the full vertical extent of 
the bed, movies of multiple independent drops, imaged 
at different, slightly overlapping vertical bed positions, 
were stitched together using the impacting sphere to align 
them horizontally and synchronize them. 

The measured intensity, /, is a function of the product 
p(pl. Here p is the density of the grain material and / the 
x-ray path length through the bed, determined from the 
cylindrical geometry of the set-up. To correct for spatial 
variations in beam intensity and camera sensitivity, cal- 
ibration curves relating / to the packing fraction cf) were 
calculated for each of the 780 x 300 pixels in the fleld 
of view [23]. Across each frame in a single movie, initial 
packing fraction varied by about 1%, indicating a uni- 
form bed packing prior to impact. From drop to drop (j)o 
varied between 0.49 and 0.52. This variation was present 
at atmospheric pressure, where the pump was discon- 
nected, as well as at reduced pressure, indicating that it 
was due to small, unavoidable differences in bed settling 
after fluidization, but not due to the evacuation of the 
chamber. 

X-ray images of the interior reveal a striking air pres- 
sure dependence of both bed and sphere dynamics (Fig. 
[1]). At atmospheric pressure (P =101 kPa), the sphere 
easily penetrates the bed, reaching the bottom of the 
system and opening up a large cylindrical hole. This 
hole closes flrst at some intermediate depth, resulting in 
a crater near the top and a gas-fllled cavity behind the 
sphere. Both of these region then flll in from the sides 
due to gravitational pressure. This process has previ- 
ously been identifled as driving jet formation [l5|, llq . 
but the role of the interstitial gas has remained contro- 
versial [17|. At P =12 kPa the overall features are still 
similar to those at atmospheric conditions, but the dy- 
namics begin to change. The cavity is smaller, and while 
the top surface of the bed still rises to compensate for the 
change in bed volume, it does not rise quite as high (Fig. 
[T^-c). Further reducing the pressure to 0.7 kPa signifl- 
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FIG. 2: (Color Online), (a) Sphere depth Zs versus time after 
impact (t = s). (b) Velocity Vs(i) computed from curves 
in (a), (c) Velocity Vs{^Zs) versus depth. All panels top to 
bottom: P = 0.7 kPa, 4.9 kPa, 8.7 kPa, 12 kPa and 101 kPa. 



cantly changes the dynamics. The bed barely rises, and 
there is signiflcant compaction in front of the sphere, ev- 
ident in the darker region under the sphere in Fig. [T]i-f. 
As a consequence, the sphere is able to penetrate only 
a short distance and the hole closes directly above the 
sphere without forming a separate cavity. 

To quantify this change in impact dynamics we plot in 
Fig. [2^ the position of the bottom of the sphere, Zs(t\ as 
it moves through the bed. At P =101 kPa the sphere hits 
the bottom with suflicient momentum to bounce back up 
a bit; at P =8.7 kPa and below the sphere is stopped well 
before reaching the bottom. From (t) we compute Vs = 
dzg/dt (Fig. [2] b, c). With deceasing ambient pressure 
there is a monotonic decrease in penetration depth and 
an increase in bed resistance. 

Inside the bed the net force on the impacting projec- 
tile is the sum of its weight, —mg^ and a drag force, Fd, 
representing the bed resistance. Models for this resis- 
tance are of the form Fd = Fc -\- c\vs\^ , where Fc repre- 
sents Coulomb friction and c characterizes the strength 
of the velocity-dependent drag. Speciflc forms include 
Fc const with /3 = 1 [5], Fc = k.\zs\ with c = [Sl,!!^, 
Fc = i^\zs\ with = 2 [^] as well as = 2 but a more 
complex Zg-dependence for Fc [lH. Note that all of these 
models predict a non-zero deceleration, a = —g -\- Fd/m, 
of the projectile and thus a Zg-dependence of its velocity 
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over the whole range from impact to final stop. 

The data in Fig. [2] for intermediate pressures reveals 
a feature not captured in these models: a region of near- 
constant velocity beginning roughly 30 ms after the im- 
pact (dotted vertical line). This behavior can be seen 
most clearly when the velocity is plotted as function 
of depth, Zg, below the free surface (Fig. [Sb). A sec- 
ond important feature is the rapid deceleration after the 
constant velocity regime, seen clearly in the traces for 
P = 4.9kPa and 8.7kPa. It is nearly as abrupt as when 
the sphere hits the bottom of the container (see traces 
for higher pressures), but occurs here sufficiently far in- 
side the bed for boundary effects to be irrelevant [lH . As 
we show below, these characteristic features are closely 
linked to the interplay between penetrating sphere, bed 
particles, and interstitial gas. 

X-ray radiography allows us to examine this interplay 
locally. In Fig. [3l we plot the change in local packing 
fraction, = (j) — (j)o^ measured along the center line of 
the path of the sphere, at three different depths, Zm below 
the surface. At P = 0.7 kPa there is a clear jump in 
well before the sphere arrives. With increasing depth this 
jump occurs further ahead of the sphere, demonstrating 
that the compacted region grows as the sphere plows into 
grains that do not flow out of the way. The situation 
is quite different at P = 101 kPa. The packing fraction 
remains constant except for a slow upturn once the sphere 
comes within about half its diameter of Zm- The width of 
this small compaction front varies little with depth. The 
overall magnitude of compaction decreases smoothly with 
increasing pressure (Fig. Hh). 

A global measure of the effect of interstitial gas is the 
rise of the top surface as the sphere burrows into the bed. 
This rise can be seen in Figs. [T^-c but is considerably re- 
duced when the system is evacuated (Figs. [Ui-f). In Fig. 

we track the level change, Sh, of the top surface for 
different pressures. At impact, the bed rapidly rises, then 
at time tm (highlighted for 101 kPa, 8.7 kPa and 0.7 kPa 
by arrows in Fig. [4^) levels off into a broad maximum of 
height Shmax and eventually falls to a final level 6h < 
(at 101 kPa this final settling does not occur until times 
much later than shown). The change in height Shmax 
increases with pressure, suggesting that the interstitial 
gas makes the bed as a whole less compressible and more 
fluid- like. 

To check this, we use x-ray images as in Fig. [T]to track 
the shape of the cavity and estimate its volume, assuming 
cylindrical symmetry. For atmospheric pressure we find 
that Sh{t) corresponds, within experimental uncertain- 
ties, to what would be expected from an incompressible 
fluid. With decreasing P the level Shmax drops below 
the value obtained from the cavity volume (double-sided 
arrows in Fig. UK), indicating a less elastic response. 

We can characterize the elasticity of the bed by com- 
paring the gravitational energy gained by the bed AUb 
to the kinetic energy lost by the sphere AKs in time 
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FIG. 3: (Color Online). Compaction front preceding sphere. 
Change in packing fraction Acj) measured along the path of 
the sphere at depths (left to right) Zm = 1-0 cm, 2.0 cm and 
3.0 cm plotted against distance from sphere tip Zs to Zm- 



tm- While at P = 0.7 kPa only 0.1% of the impacting 
sphere's energy is transferred to the bed, this value in- 
creases to about 45% at atmospheric pressure (Fig. 
This can be compared to values around 10% found in 2D 
simulations without interstitial gas [l2[. 

To understand the fluid-like behavior at large P, we 
examine the rate of gas flow through the bed. If the 
timescale for the expulsion of the gas from the bed is 
significantly longer than the timescale for the granular 
flow, then gas trapped and compressed by the bed can 
create pressure differences capable of supporting the bed 
|ld llSl, l20|, [2l|. From Darcy's law and the continuity 
equation for the gas flow, one can derive a diffusion equa- 
tion for the gas pressure ^ = P^§^, with diffusion con- 
stant D = ^^Y^-rAV where /i is gas viscosity and k bed 
permeability pj|, |20|. For our experimental conditions 
D ^ 5 cm^/s at P = 101 kPa. The timescale to dif- 
fuse across the depth of our bed (a distance L = 8.5 cm) 
is td = / D ~ 140 ms, significantly longer than the 
time tm ~ 30 ms for the bed to rise to Shmax- This 
suggests that air trapped in the bed interstices prevents 
compaction at large P. Since the permeability depends 
on the grain diameter according to k ^ (P^ this cushion- 
ing effect would be less pronounced with larger grains, 
as noted in [3]. Conversely, we expect the behavior of 
larger grains to resemble that found at our lowest pres- 
sures. Indeed, the trajectory of the sphere at 0.7 kPa is 
qualitatively similar to trajectories measured by Durian 
et al. in 250/im - 350/im glass spheres [6,,^]. 

The air-mediated response of the granular bed directly 
affects the drag on the impacting sphere. Comparing 
Figs. [4^ and [2)3,c we see that the onset of the constant 
velocity regime coincides with tm^ the start of the plateau 
in bed level rise. During this stage the material displaced 
by the sphere is flowing mostly into the cavity behind it. 
Consequently, the sphere is not affected by the full bed, 
but instead by a more local region. For Vs ^ 1 m/s, 
as for our data, the drag (ppD^vl is within a factor of 
two of the weight of the sphere. This suggests a large 
reduction of Fc in the force law in the observed constant 
velocity regime. When the top level of the bed begins to 
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FIG. 4: Bed dynamics, (a) Rise of bed surface 6h. Arrows 
mark tm for three pressures. Double-sided arrows denote rise 
needed to conserve bed volume . Resolution 6h was limited 
by pixel size to ~0.5mm. (b) Maximum change in packing 
fraction in front of sphere (Q) and ratio of potential energy 
needed to raise bed by dhmax to kinetic energy lost by sphere 
at time tm (■)• Error bars for A(/)max are due to fluctuations 
in A0, and for AUb/^Ks due to uncertainty in tm- 

fall again, the cavity behind the sphere has pinched shut 
(Fig. [lb), trapping an air pocket below the surface [l6|. 
With the falling bed and trapped air pocket, the bed 
material is no longer free to flow out of the way of the 
sphere. Coulomb friction is again set by the full weight 
of the bed, resulting in an increase in Fc that quickly 
decelerates the sphere and brings it to rest. 

As a result, a friction term of the form Fc = k,\z\ can- 
not capture the full range of observed behavior, even if 
Hi is made to depend on pressure. Such pressure depen- 
dence was very recently proposed by Cabarello et al. 
who measured the trajectory of a sphere impacting 40 fim 
sand by tracking a string attached to it. Our results agree 
with their conclusion that the shallower penetration at 
lower pressures is due to increased friction. However, 
since A(j) > (Fig. [3]) at all pressures, drag reduction 
is not simply due to fluidization of grains ahead of the 
sphere. Instead, our measurements of the local packing 
fraction indicate that in the presence of air the bed as a 
whole behaves more like an incompressible fluid, allowing 
the impacting sphere to penetrate deep and create a large 
cavity. In the absence of air the bed compacts much more 
strongly ahead of the sphere, rapidly dissipating energy 
and decelerating the descent. 
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